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NA44 uses a 512 channel Si pad array covering 1.5 < n < 3.3 to study charged 
hadron production in Pb+Pb collisions at the CERN SPS. We apply a multircso- 
lution analysis, based on a Discrete Wavelet Transformation, to probe the texture 
of particle distributions event-by-event, by simultaneous localization of features in 
space and scale. Scanning a broad range of multiplicities, we look for a possible 
critical behaviour in the power spectra of local density fluctuations. The data are 
compared with detailed simulations of detector response, using heavy ion event 
generators, and with a reference sample created via event mixing. 



An ultrarelativistic collision of heavy ions presents a phenomenon whose 
most interesting features are conditioned by the large multitude of degrees of 
freedom involved, and yet offer an opportunity for the fundamental physics of 
the strong interaction to manifest itself. The very notion of a phase transition 
in such collisions is inherently of a multiparticle nature. Truly multiparticle 
observables, defined on event-by-event basis - a few have been constructed so 
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far - therefore attract attention. Recently published event-by-event analyses 
of the 158 GeV/A Pb + Pb data either deal with a small number of events 111 
or analyse properties of a large ensemble oLevents using a single observable 
(px) to compare different ensemble averagesa. In the first case, accumulation 
£if feature information from large data sets remains open. In the second case 
□, one can not establish a scale independency in event textures by observing 
a logical consequence thereof 0. Furthermore, an ensemble average on a set 
of post-freeze-out events is not representative of the pre- freeze- out history of 
those events, due to the dramatic non-stationarity of the open system, with 
a consequent lack of ergodicity. 

Here we concentrate on texture, or local fluctuation observables, where a 
single event determines its own correlation/fluctuation content, and the scale 
composition thereof manifests itself in the observables in a positive way. The 
idea to look at particle_distributions in rapidity y to search for critical be- 
haviour was proposed au based upon a Ginzburg-Landau type of multihadron 
production theory □, where a random hadronic field <fi(y) plays the role of an 
order parameter in a hadronization transition. Enhanced large scale correla- 
tions of hadrons in y at the phase transition would_signal critical fluctuations 
in the order parameter. Stephanov and coworkers □ indicated a second order 
QCD phase transition point which should exist under certain initial condi- 
tions, within the reach of today's experiments. 

In our work, a power spectrum analysis of event texture in pseudorapid- 
ity rj aad azimuthal angle £, based on a Discrete Wavelet Transformation 
(DWTJ2, is performed on a number of large event ensembles sampled accord- 
ing to their multiplcity, thereby studying the impact parameter dependence of 
the observables. DWT quantifies contributions of different £ and rj scales into 
the overall event's texture, thus testing the possible large scale enhancement. 

The SPS beam was collimated to a 1 x 2 mm profile. The NA44 Si pad 
array, installed.-! cm downstream from the target, in the magnetic field of 
the first dipoleQ, measured ionization energy loss of charged particles in its 
512 300 /im thick Si pads. The silicon detector had inner radius 7.2 mm and 
outer radius 43 mm. The detector was split radially into 16 rings of equal rj 
coverage. Each ring was divided azimuthally into 32 sectors of equal angular 
coverage to form pads. The pads were read out by AMPLEX □ chips, one 
chip per sector. (S-electrons, produced by the Pb beam traversing the target, 
were swept away to one side by the dipole magnetic field (< 1.6 Tl). Only 
the (5-electron-free side was used in this analysis. Channel pedestals had, on 
the average, FWHM = 0.48 < dE > of 1 MIP. In the texture analysis, every 
event was represented by a 2D array of the calibrated digitized amplitudes of 
the channels ( an amplitude array). Empty target runs were used to measure 
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the background. Cross-talk in the detector was evaluated off-line. 

DWT formalizes the images of the PbPb collision events in pseudorapidity 
r] and azimuthal angle £ by expanding them into a set of functions orthogonal 
with respect to scale and position, and allows one to accumulate the texture 
information by averaging the power spectra of many events. The simplest 
DWT basis is the Haar one, built upon the scaling function <f>{x) = 1 for 
< x < 1 and otherwise. If the interaction vertex lies on the detector's 
symmetry axis, every pad's acceptance is a rectangle in the (£, rf) space. Then, 
the Haar basis is the natural one, as its scaling function in two dimenstions 
(2D) <!>(£, 77) — 4>{Q(j){ri) is just a pad's acceptance (modulo units). We set up 
a 2D wavelet basis: 

*i., JC V) = 2 m * A (2 m C - i, 2 m V - J) (1) 

&m,i,j (Ci V) i s constructed from $(£,77) similarly. Here, m is an integer scale 
fineness index; i and j index the discrete positions of pad centers in £ and 77 
( 1 < to < 4 and 1 < i,j < 16 because we use 16 = 2 4 rings and 16 sectors 
). Different values of A (denoted as £, 77, and £77) distinguish, respectively, 
functions with azimuthal, pseudorapidity, and diagonal texture sensitivity: 

# C = ^(C)M), = #0tK»7), ^ = HCWv) (2) 

In the Haar basis, for any variable x 

ip(x) — {+1 for < x < — ; — 1 for — < x < 1; otherwise} (3) 

is the wavelet function. Then, i ■ with integer to, i, and j aie known H 
to form an othonormal basis in L 2 (R 2 ). We adopt the existing t3 ID DWT 
power spectrum analysis technique and expand it to 2D. Track density in 
an individual event is /?(£, 77) and its local fluctuation in a given event is 
a 2 = (p — p, p — p) , where p is the average p in the given event. Using 
completeness of the basis, we expand 

p-p = {p, iW )¥* ^ - (p, ^.,„, (4) 

Notice that /5, being constant within detector's rectangular acceptance, 
is orthogonal to any ^ 3 with m > 1 . Due to the orthonormality condi- 
tion (^m,i,j>^m',i',j') = ^A,A'^m,m"^i,i'^-,i' , the p -p components for different 
scales do not form cross-terms in the a 2 sum, and the sum contains no cross- 
terms between p and p for the four observable scales. Instead of a (p, <fr m =5.i,j) 
set, the amplitude array - its closest experimentally achievable approximation 
- is used as the DWT input. 
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Figure 1. Power spectra of 7 X 10 3 events in the multiplcity bin 379 < dN/ drj < 463. O 
- true events, A - mixed events, □ - the average event. 



The Fourier images of ID wavelet functions occupy a set of wave numbers 
whose characteristic broadness grows with scale fineness mas2 m ; 2 2m should 
be used in the 2D case. In 2D, we find it most informative to present the three 
modes of a power spectrum with different directions of sensitivity P^ T, (m), 
pt(m), P v (m) separately. We define power spectrum as 

^M = ^E^><m> 2 > (5) 

where the denominator gives the meaning of spectral density to the observ- 
able. So defined, the P(m) of a random white noise field (for example) is 
independent of m UJ. In the first approximation, the white noise example 
provides a base-line case for comparisons in search for non-trivial effects. 

We used WAILI software library to obtain the wavelet expansions. 
Figure [j] shows the power spectra for one multiplicity range. The first striking 
feature is that the power spectra of physical events are indeed enhanced on the 
coarse scale. The task of the analysis is to quantify and, as much as possible, 
eliminate "trivial" and experiment-specific reasons for this enhancement. 

The average event, formed by summing amplitude arrays of the measured 
events within a multiplicity range, and dividing by the number of events, has 
a much reduced texture as fluctuations cancel. However it retains the texture 
associated with the d 2 N/ drjdC,, with the dead channels and the finiteness of 
the beam's geometrical cross-section. A better way to get rid of the "trivial" 
texture is to use mixed events. The event mixing is done by taking different 
channels from different events. Therefore, the mixed events preserve the tex- 
ture associated with the detector position offset, the inherent dN/ drj shape 
and the dead channels. In order to reproduce the electronics cross-talk ef- 
fects in the mixed event sample, mixing is done sector-wise, i.e. the sectors 
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Figure 2. Multiplicity dependence of the texture correlation. — the NA44 data, • - 
RQMD. The boxes show the systematic errors vertically and the boundaries of the multipl- 
city bins horizontally; the statistical errors are indicated by the vertical bars on the points. 
The rows correspond to the scale fineness m, the columns - to the directional mode of the 
power spectrum A (which can be diagonal £77, azimuthal f, and pseudorapidity rf). 



constitute the subevents subjected to the event number scrambling. In other 
words, the mixed events preserve the texture coupled with the channels of 
the detector. This is static texture as it reproduces its pattern event after 
event; we are interested in dynamic texture. We reduce sources of static 
texture in the power spectra by empty target subtraction and by subtrac- 
tion of the mixed events power spectra, thus obtaining the texture correlation 
P x (m)true — P X {fn)mixed- Its multiplicity dependence is plotted on Figure 
|. For comparison with models, a MC simulation (done with RQMD E3) in- 
cludes the known static texture effects and undergoes the same elimination 
procedure. This allows the effects irreducible by the subtraction methods to 
be taken into account in the comparison. One such example is the finite beam 
size, which has been shown by the MC studies to cause the RQMD points to 
rise with dN/ dr\. 
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The systematic errors were evaluated by removing the Pb target and 
switching magnetic field polarity to expose the analyzed side of the detector 
to (5-electrons, while minimizing nuclear interactions. All correlations (i.e. 
deviations of P x {m)t rue from P x {m) m ixed) in such events are considered to 
be systematic errors. Thus this component of the systematic error gets a sign, 
and the systematic errors are asymmetric. The other component (significant 
only on the coarsest scale) is the uncertainty of our knowledge of the beam's 
geometrical cross-section. 

This novel method of event-by-event analysis, applied to the SPS PbPb 
data, does not reveal any evidence of critical phenomena. 
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